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they differ only in the degeneracy of their eigenvalues. This is true for finite 
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the higher dimensional ones. The Bethe states are highest weight states of 
the quantum group, except for some states with energy zero. 
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1 Introduction 



The recent interplay between the field of solvable two-dimensional lattice ( or quan- 
tum spin chain ) models and quantum groups, has generated a lot of interesting 
results. One particular way of building models, which are quantum group invariant, 
uses the Temperley-Lieb (TL) algebra Q, satisfied by the Hamitonian density U k : 

Ul = yjQU k , U k U k±1 U k = U k , 

[U k ,Ui] = 0, \k-l\> 2. (1.1) 
The Hamiltonian is now given by the following sum over N sites: 

N 

#(?) = £ £4. (1.2) 

k=l 

The paper is organized as follows. In section 2, we describe the representations 
of the TL algebra, constructed as projectors on total spin zero of two neighbouring 
spins. In section 3, we discuss a modified coordinate BA and show the reasons, why 
the techniques develloped for the spin 1/2 XXZ model don't work here. Section 3 
contains the algebraic details of the computation and section 4 is reserved for the 
conclusions. 

2 Representations of the Temperley-Lieb algebra 
as spin zero projectors. 

Representations of the TL algebra, commuting with quantum groups, can be con- 
structed in the following way|Q. Suppose U q {X n ) is the universal envelopping alge- 
bra of a finite dimensional Lie algebra X n , equipped with the coproduct A : lA q — > 
U. q ®U q || . If now 7r :U q —* EndV^ is a finite dimensional irreducible representation 
with highest weight A and we assume that the decomposition Va <S> Va is multiplic- 
ity free and includes one trivial representation on Vq, then the projector Vq from 
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V\ ® V\ onto Vq is a representation of the TL algebra. The deformation parameter 
q, which plays the role of a coupling constant in the Hamiltonian, is related to Q 
as: 



where p is half the sum of the positive roots. 

By construction Vq commutes with the quantum group U q (X n ). 

Since we are not going to use any group-theoretical machinery, we will just lift 
the relevant formulas off Batchelor and Kuniba[|J in order to display explicitly the 
Hamiltonians to be diagonalized. 

We will consider the following specific cases, (V\,U q (X n )) = (V2skiMq{Ai)) for 
spin s, (Vad U q (B n ),(V'A. 1 , U q {C n ) and (V^, U q (D n ). I.e, we treat the g-deformations 
of the spin-s representation of sl(2) and the vector representations of so(2n + 
1), sp{2n) and so{2n). V\ denotes the U q (X n ) module with highest weight A. Ai is 
a highest weight of X n . 

Introduce the following notation. Let e^, i — 1, . . . , n be orthonomal vectors and 
express the fundamental weight, II = Aj + . . . + A n , the set A of weights and the 
coupling constant ^[Q = — 2A as: 




(2.1) 



A l : A 



{s(e x - e 2 ), (s - l)(ei - e 2 ), . . . , -s(e x - e 2 )} 



Ai = (ex - e 2 )/2, 



p=(e 1 - e 2 )/2, 



J = {s,s 



-1,...,- 



s}, 



yjQ= [2s + l]; 
B n (n>2): A= {0, ±e 1; . . . , ±e n }, 



A« = t\ + . . . + e i: (1 <i < n), 
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= (ei + ... + e„)/2, (i = n), 
p=(n- l/2)ei + . . . + e n /2, 
J={0,±l,...,±n}, 

= (-1) A , 
^/Q = [2n-l][n + l/2]/[n-l/2]; 

C n :.4 = {±ei,...,±e n }, (2.2) 

Aj = ei + . . . + ej, (1 < i < n), 
p = nei + . . . + e n , 
J={±l,...,±n}, 
e(p) = sign(p), 
^Q = [ra][2ra + 2]/[n + l]; 
D n :A = {±ei,...,±e n }, 

Ai = e 1 + . . . + d, (1 < i < n — 1), 
= (ei + . . . + e n _i - e n )/2, (i = ra-l), 
= (e 1 + . . . + e„_i + e n )/2, (i = n), 
p=(n- l)ei + . . . + e„_i, 
J = {0,±l,...,±n}, 

= [2ra-2][ra]/[ra-l]; 
e(p) = 1. 

For p e J the symbol /2 is defined as /i = + (1 ± l)/4 for A 1 with s e Z + (1 ± l)/4 
and /i = with the exception of = 1 for B n . The g-number notation is [x] = 
(q x — q~ x )/ (q—q^ 1 ) ■ For X„ = B n , C n , D n , we extend the suffix of e M to — n < p < n 
by setting e_ M = — e M (hence eo = 0). Using the index set J, above, we can write 
A = {p{ei - e 2 )} for A x and A = {e^p E J)} for B n , C n , D n . 

Denoting by E^ u e EndV\ the matrix unit, having all elements zero, except at 
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row n and column u, the projector can be written as 

*Po = Q- 1/2 E <n)<v)q- <e » +e - p> E, u ® E^ v . (2.3) 

In the following we will refer to all models generically as higher spin models for 
simplicity, even when not talking about A\ 

Consider a one-dimensional chain of lenght N with a "spin" at each site. The 
spin variables range over the set of weight vectors v^/i G J and our Hilbert space 
is an iV-fold tensor product V\ <g) . . . <g) V\. For A\, these are the g-analogs of the 
usual spin states. 

The Hamiltonian densities acting on two neighboring sites are then given by: 

(k,l\U\i,j) = e(i)e(k)q-< e ^>5 l+3fi 5 k+lfi . (2.4) 

3 The coordinate Bethe - Ansatz 

All the above Hamiltonians are U(l) invariant and we can classify their spectra 
according to sectors. For A\(s = 1) the commuting operator is the total spin 
S z = J2k=i $k an d we set the conserved quantum number r = N — S z . In general 
it equals r = N * uo — S z for A\ and B n and r = iV * (uo — 1/2) — S z for C n and D n . 
We set uj = max J. 

Therefore, there exists a reference state |fi), satisfying H \Q) = E \£l), with 
E = 0. We take \Q) to be \Q) = Uk |w, k). 

In every sector r there are eigenstates degenerate with \Q). They contain a 
set of impurities. We call impurity any state obtained by lowering some of the 
|o;,A;)'s, such that the sum of any two neighboring spins is non-zero. Since H(q) 
is a projector on spin zero, all these states are annihilated by H(q). In particular, 
they do not move under the action of H(q), which is the reason for their name. 

We will now start to diagonalize H(q) in every sector. Nothing interesting 
happens in sector r < 2u. Sector r = 2u is more interesting, although still trivial, 
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since it contains one free pseudoparticle. The main result of this paper is to show- 
that H(q) can be diagonalized in a convenient basis, constructed from products of 
single pseudoparticle wavef unctions. The energy eigenvalues will be parametrized 
as a sum of single pseudoparticle contributions. 



3.1 The sector r = 2cj, containing one pseudoparticle 

Starting with r = 2u>, we encounter the situation, where the states \j, k) and 
|— j, fc±l), j 7^ uj occur in neighboring pairs. They do move under the action 
of H(q) and mix with states containing one \—u,k). Eigenstates are a superpo- 
sition of |x[-a;]) = ( . . . ujujojuj —ujujuj . . .) and \x[j,~j]) = (. . . uuu +j — jujuj . . .), 



\2u; ...) = J2i a M \ x l~^) + 12 h A x ) \ x U-j])}, (3.1) 
x j 

where Y!j means j E J* = J — {±c<j} and the ellipses stand for parameters the 
eigenvector is going to depend on. When H(q) now acts on \2u>; . . .) it sees the 
reference configuration, except in the vicinity of x and we obtain the eigenvalue 
equations 



(E - g 2 <^> - g - 2 <w>) au ( x ) = au ( x + !) + a ^ x _ !) + 

J2e(u;)e(l)q- <e « +ei ' p> b t (x - 1) + ]T e(-u)e(l)q- <e - u+ei ' p> b t (x) 
i i 

Ebj(x) = e{uj)e{l)q- <e > +e "' p> auJ {x + 1) + e{-u)e{l)q- <e ^ e -' p> a u {x) + 

j2<Mi)q- <e3+euP> bi{x), jer. (3.2) 
i 

Eliminating the b/s, we get an equation very similar to the XXZ model: 

q~ 2<ei ' P> )aM = a u (x + 1) + a w {x - 1). (3.3) 
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We will treat periodic boundary conditions maintaining translational invariance 
in the following sections. They demand a w (x + N) = a w (x) and bj (x + N) = bj(x). 
We parametrize as: a u (x) = a^ x and bi{x) = b£ x , I G J*. Substituting this into 
equ.( p.2| ) we get two eigenstates and their energies 

a w = e{-uj)q <e ^ p> + e{u)q- <e ^ p> C l =T{C l ) (3-4) 
ft, = e{l)q- <ei " p> , / G J* (3.5) 

E x = E^ <£i ' p> + r (Or(r 1 ) = E^ 2<ei,p> + e + r 1 

l lej 
and a highly degenerate solution with E2 = 0, with the following constraint on the 

parameters: 

E eilX^q-^+^k + r(Oe(^)g- <e -^ > a L , = 0. (3.6) 
1 

Here ^ = e l9 , 8 being the momentum determined from the periodic condition to be: 
9 = 2nl/N, with I integer. 

We describe this situation by saying that we have two types of pseudoparticles 
with energies E 1 and E 2 . Whereas the pseudoparticle \2u;9) 2 is degenerate with 
i.e. propagates with energy E 2 = 0, the pseudoparticle \2u;9) 1 , propagates 
with energy 

E x =-2A + 2cos#, 2A = -J2q' 2<ei ' P> - (3.7) 

As mentioned before, the energy eigenvalues are going to be parametrized as a sum 
of single pseudoparticle energies. Thus we write: 

e = I>„d> <ei ' p> + r (£n)r(C 1 )), (3.8) 

n=l I 

where e n depends on which pseudoparticle we use: e n = 1 for E = Ex and E = 
E 2 = 0. 

3.2 Two pseudoparticles and the XXZ Bethe - Ansatz 

The next higher sector would be r = 2u + 1, but let us treat r — 4u first, since 
then we can compare it with the first nontrivial sector in the XXZ model. 
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This sector contains states, which consist of two interacting pseudoparticles. We 
seek these eigenstates in the form: 



|4w; • • -) eie2 = E {^uu(x 1 ,x 2 ) \x 1 [-u],x 2 [-u]) +J2b U}i (x 1 ,x 2 ) \x x [-uj},x 2 \i-i}) + 

xi<X2 i 

^2b juJ (x 1 ,x 2 ) \xi[j,~j],x 2 [-w]) + ^2^2b ij (x 1 ,x 2 ) \x 1 [i~i],x 2 [j,-j])}. (3.9) 

3 » 3 

Translational invariance now specifies a^i^i, x 2 ) = C, Xl a UJU j(n) and similarly for the 
other wave functions, where n = x 2 — X\. Periodic boundary conditions require that 

b tuJ (n) = Fbrt(N-n),b ij (n) = ?b ji (N-n), (3.10) 

where £ = £i£ 2 (£» = i — 1,2) and the total momentum is 6i + 6 2 = 2nl/N, 
with I integer. 

According to equ.( |3.8|) , we will parametrize the energy as 

e= E^Eff-^ + rcwrc^ 1 )]. (3.11) 

n=l / 

Let us take the block e\ = e 2 = 1 first. We try to build 2-pseudoparticle 
eigenstates out of translationally invariant products of 1-pseudoparticle excitations 
at x\ and x 2 with weight functions Di(xi,x 2 ), i — 1,2: 



\4u;9 u e 2 } 11 = E { J Di(x 1 ,a; 2 )[r(£r 1 )|x 1 h-])+E e (0^ <ei ' P> ki[^])] 

[ rfe- 1 )]^^]) + E<ik~ <e ^ > i^w]>] + 

D2(xi > x 2 )[r(e 2 - 1 )|x 1 [- w ]> + Ee(;)?" <ei,p> l^ib-,-i]>] 

j 

[rcer 1 )!^!— 1> + E e W^ <e! ' p> k 2 [^])] }• 
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Comparing this with equ.( |3.9| ) and using translational invariance, implying ^(n) 
CD^N -n), we get 



b iuJ ( n ) 
b i:j (n) 



c(<) q-^im 1 ) D 1 (n) + rfe" 1 ) D 2 (n)}, 
e^q-^im^D^n) + )D 2 (n) ], 

e(i)eO') g- <6l+ej ' p> D(n), 3 < n < JV - 3, 



(3.12) 



where D(n) = Di(n) + D 2 (n). 

Applying -ff(g) to the state of (|3.9|), we obtain a set of coupled equations for 
a>uxj(n), hj {n). Following 0, we split the equations into far equations, when excita- 
tions do not meet and near equations, containing terms when they are neighbors. 
The far equations are: 

(E - 2q~ 2<e ^ > - 2q 2<e ^ > )a ulU! (n) = (1 + O' 1 )^ + !) + (! + £K> - 1) 



J2 e{l)e{u)q-<^ +e ^>[C l bUn + 1) + b ul (n - 1) ] + 
l 

^eilH-^q-^-^ibiM+b^n)}, 2<n<N-2, 
i 

(£ _ g -2<e„,p> _ ? <2 W >)^.( n ) = + 1} + ^.( n _ 1} + 

e(j)e(-o;)g- <e ^ e -' ,> o^n) + eO>(u;)<r <e ' +<w> a wt > + 1) + 



(3.13) 



<ei,p> 



e{uj)q- <e «> p> C%{n + 1) + e(-cu)g <e ^ > ^(n) + 



2 < ra < iV- 2, 



(3.14) 



( E - g~ 2 <^> - q 2 <^>)b, JUJ {n) = b juJ (n - 1) + b uj (n + 1) + 
e(j)e(-a;)g- <e ^ e - p> a^(n) + eO>(uO<T <e ' +<w> £a«*,(ri - 1) + 



<ei,p> 



e(uj)q-<^>b 3l (n - 1) + e^g^^n) + e(j)<T <ej ' P> ^ W 



/ 

3 <n < N -3, 

Ebijin) = e(^)g- <e -" > [e(7)g- <e -^ > ^(n - 1) + e{j)q- <e » p> b M {n + 1)] + 
e (_ w ) g <^P> [eWg -<^>^. (n) + e (j)g-< e ^>^(n)] + 



(3.15) 
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/ 



J2<^~ <ei ' P> Wh~ <ei,P> bij(n) + eU)q- <e ^ p> b a {n)], 
i 

3<n<iV-3. (3.16) 

We already know them to be satisfied, if we parametrize Z?i(n) and D 2 {n) by 
plane waves: 

Di(n) = D 2 {n)=$Q. (3.17) 

The real problem arises of course, when pseudoparticles are neighbors, so that 
they interact and we have no guarantee that the total energy is a sum of single 
pseudoparticle energies. 

Let us now have a lightning review of the XXZ coordinate BA in order to be 
able to comment on the features, which are not going to survive generalizations to 
the present models. 

The equations in the sector r X xz = 2 are 

(E - 2q^ - 2g - 2) a(x x , x 2 ) = (3.18) 

a(xi + 1, x 2 ) + a{x\ - 1, x 2 ) + a(x 1 , x 2 + 1) + a(xi,x 2 - 1), (3.19) 
if x\ and x 2 are not neighbors. In case they are we get 

(E — — q~^) a(xi, X\ + 1) = a{x\ — 1, x 1 + 1) + a(xi, Xi + 2). (3.20) 



One now supposes, that the parametrization equ.( |3.9| ) for s = 1/2 solves both the 
above equations. In this case we are allowed to set x 2 = x\ + 1 in equ.( |3.18D and 
subtract it from equ. (|3.20| ), yielding the following consistency condition : 



- {q 2 + q 2 )a(xi,xi + 1) = a(xi,xi) + a(xi + l,X\ + 1). (3.21) 

This gives the BA equation for the XXZ model, determining the 2-body phase 
shift: 

= + + (3.22) 
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This type of procedure only works for the spin 1/2 XXZ model, due to the 
following fact. When the two pseudoparticles come together in a configuration like 

(. . . + + H h + + ...) and when H is applied to the two down spins, it gives zero, 

since their total S z equals —1. But whenever two excitations approach each other 
becoming neighbors and the Hamiltonian applied to them yields a nonvanishing 
resul t, then the representation like equ.( 3.12j ) cannot solve both the far and near 
equations. Yet this is exactly the situation arising for higher spins. As we shall 
see, in this case, the representation equ.( p.l^ ) has to modified |TTJ], the two-body 
wavefunction developing a "discontinuity" at minimum separation. We call this the 
spin zero rule. 

Now back to our problem to solve the near equations. They are: 



(E-q 



-2<e„,p> 



g 2<e "' p> )a^(i) = (i + _1 ^(2) + 



J2 t{l)q- <e ^ p> [t{uj)q' <e ^ p> r'M2) 



(E - q-^^b^l) = rXj^) + 
e{j)e{-u)q- <e ^> 0^(1) + e(j)e(u)q~ <e ^ p> 0ww (2) + 



<ei,p> 



-<e w ,p> t~ 1 



r\(2) + e(j)q-< e ^>b^(l) 



(£ - g 2<e -' p> )M 2 ) = M3) + 

e(j)e(-u)q-<^ p> aU^) +e(Mu)q- < ^ p> ^aU^ + 



<ei,P> 



,(_ w ) g <w> 6 . i(2 ) + e ( 3 ) q -<^P>b lu) (2) 



EbM) = e(u)q 



-<e„,p> 



e(i)q 



-<ei,p> 



^ j (l) + e(j)q-< e ^>bU^} + 



e{-u)q 



<e^,p> 



i 



-<e;,p>[ 



[e(i)g-< B *.">6 wi (2) +e(j)q- <e ^ > b iw (2)\ + 
e^q-^b^+e^q-^k^)], 
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Eb u (2) = e^e^q-^+^^l) + 6^(3)] + 

e(i)e(-a;)g- <ei - e -' )> [6 1Ji (2) + b iu (2)] + ^ e(/)g- <e ^'" > [^(2) + b u (2)] + 

i 

e(-i)e(Og" <e!_ei,p> Bf ) . (3.23) 

i 

Here some new states are showing up. Bf are the wavefunctions of the states 
of the type (. . . toui I — I ojuj ...),/ ^ i. Applying H(q) to them we obtain the 
system: 



( £7 _ 9 2<e„p> )jB (0 = J- e(/) e (j>< e ' +e -^Sf ) + e (/) e (-0g- <ei - e -' 5> ^(2), (3.24) 
yielding 



6 (o = ,j i ^ < ' , :::, M ^ 



(E-Y. K j, j ?-,<l- 2< '"* > 
Eliminating B^P from equation ( |3.23|) , we get 

E + 2A + + q ^lp> U2) = <^<^ <ei+ ^ P> \^M) + M3)] + 
e(z)e(-^)g- <e - e -" > ^(2) +6 to (2)] + 

j2<lXi)q- <e * +e " P> ibu(2) +6«(2)]. (3.26) 

In order to solve these equations, we follow [jDj] and now leave the value of the 
wavefunctions for nearest separation as arbitrary parameters: 

Mi) = ^W^ <ei ' p> [r(r 1 ) J Di(i) + r(r 1 ) J D 2 (i)] + ^(i), 

«U(2) = e «g-< e -^>[r(r 1 ) J Di(2) + r(r 1 )£> 2 (2)] +^„(2), 

6y(2) = e(z)e(j)q-<^>D(2) + ^(2). (3.27) 

In order for this modification to leave the far equations still satisfied, the fol- 
lowing conditions have to hold: 

(1 + r 1 )^!) + E 6(/)eMg-< e ' + -^>^(l) + 
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Y,e(l)e(-uj) q -<^-^>F b j2) = 0, 



^ buj (l)+J2e(l)e(-u)q-< e ^>^ bl .(2) = 0, 
i 

^(2) + y te(l)e(u)q-«* +e «><»,F b . l (2) = 0. 



(3.28) 



Now using equs.( |3~l~2D and equ. (|3.27|) in equ. (|3.26|) , we get the following equation 
for ^(2): 

(3.29) 



^(2) = ^^r, ie-l ■ 



E + 2A 

Doing the same with equ. (|3.23|) , we get: 



F btJ (2) = 0, i^j. 



(3.30) 



These results for jF fe ..(2), jF feij (2) are reasonable. In the first case the colliding exci- 
tations satisfy the zero spin rule and we get a non-zero result, wheras in the second 
case the rule is not satisfied and we get zero. Using this in equ.( [3.28] ), we obtain 
for the remaining constants: 



-2<e LJ ,p> _|_ q2<e LJ ,p> 



+ 2A)F bii (2), 



^„(1) 

F K .{1) = -e{j)e{-uj)q-< e ^> T bii {2\ 
^(2) = —e(j)e(u)q~ <e ' +Bu ' p> T bii {2) , j e J * 



(3.31) 



Substituting finally the complete parametrization into the remaining near equatins, 
get the following Bethe-Ansatz equation: 

D(2) £ 



E + 2A 1 + f 



D(l), 



which can also be rewritten as 



N 



e 2 [(i+r 1 )6-2A- J B] 



(3.32) 



(3.33) 



52 ^ ei[(i + e~ 1 )6-2A- J E;]- 

Using the explicit form of the energy, the set of equations determining the spectrum 



are: 



N = l + & 2 -2A6 N = 



(3.34) 
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Notice that this equation is independent of n or any other representation specific 
quantities. All the models considered show therefore an equivalent spectrum, when 
parametrized in terms of A. 

In particular, this is the same consistency condition one finds for the XXZ model, 
showing that for e\ — e 2 — 1, even for periodic boundary conditions, the spectra 
of all our models are equivalent to the spectrum of the XXZ model, if expressed in 
terms of A . 

We will refrain from discussing the other two blocks: £1 = 1, e 2 = and t\ = 
e 2 = 0, since the calculations are analogous to the ones presented above. For details 
the reader might consult reference ||1 1|| . Anyhow, the first of the two cases doesn't 
show up for the more interesting situation of free boundary conditions. It is too 
asymetric to satisfy free boundary conditions. This is the reason, why for free 
boundary conditions, the spectrum of all of our models is equivalent to the XXZ 
spectrum. The block ei = e 2 = has E = and the BA equation reduces to £ N = 1, 
being highly degenerate. The eigenvalue E = also occurs in the XXZ spectrum, 
albeit with different degeneracy. 

3.3 One pseudoparticle and impurities 

Since the setup with pseudoparticles and impurities is a little different from the 
case of two pseudoparticles, we will dedicate some space to it. 
The eigenstates sought for, will be like: 



2u + 1;...) = { a ^k(x 1 ,x 2 ) \xi[-u], X2[k\) + ^2b jk (x 1 ,x 2 ) \x 1 \j,-j),X 2 [k]) + 



Xl<X2 




(3.35) 



j 



Translational invariance and periodic boundary conditions impose: 




(3.36) 
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a„ k (n) = CaUN - n), b jk {n) = Chj(N - n), (3.37) 

where n — x 2 — Xi, £ = £1^2 (6 = e i£, % i = 1,2) and the total momentum is 
1 + 2 = 2nl/N, with I integer. 

Let us take the block e± — 1 , building eigenstates out of translationally in- 
variant products of 1-pseudoparticle excitations at x\ and an impurity at x 2 with 
weight functions D i (x 1 ,x 2 ), i = 1,2 as in the previous section. This yields the 
parametrizations: 

a„ k (n) = T(^ 1 )D 1 (n), 
aUn) = r(£ 1 )D 2 (n), 
b jk (n) = eU)q- <ei ' p> D 1 (n), 

hj(n) = e(j)q- <ei ' p> D 2 (n). (3.38) 
The far equations for the impurity at the right are now: 

(E - 2q- 2<e »* > - 2q 2<e ^ > )a U!k (x 1 ,x 2 ) = a uk ( Xl - 1, x 2 )) + a wfc (a;i + 1, x 2 ) + 

^(Oe^-^+^M^ - 1,2:2) + i2e(l)e(-u)q- <e ^ > b lk (x 1 ,x 2 ), 

1 1 
xi + 2 < x 2 < N-xi -2, (3.39) 

Eb jk ( Xl ,x 2 ) = e(j)e(uj)q^ +e -' p> a ujk (x 1 + l,x 2 ) + 

e(j)e(-w)g- <e ^"' p> a cjfc (a; 1 ,x 2 ) + 

E e(j)e(l)q- <e ^> [e(i)q- <ei ' p> b lk (x u x 2 ), 
1 

xi + 3 < x 2 < N - xi - 3, (3.40) 

and analogous equations for the impurity at the left. Eliminating the 6-functions, 
we get: 

(E-J2 Q~ 2<ei ' P> ) a u J k(x 1 ,x 2 ) = a wfc (xi -l,x 2 ) + a^xx + 1, x 2 ), 
leJ 

(E - E Q~ 2<ei ' p> )aku,(xi, x 2 ) = a klM (x 1 ,x 2 - 1) + a ku (x 1 ,x 2 + 1), 
leJ 

xi + 3 < x 2 < N-X! -3. (3.41) 
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We know them to be satisfied, if the energy is given by equ.( 3.11 ). The near 
equations require of course the by now costumary treatment of modifying the ansatz 
of the wavefunctions at nearest separations. 

The near equations for the impurity at the right are: 

(E-2q- 2<e - p> )a uk ( Xl x + l) = 
i 

a uk {x - l,x + 1)) + Y,< l )<u)l~ <ei+eu ' P> bik{x -l,x+ 1), (3.42) 

i 

Eb jk (x, x + l) = t{])t{uj)q- <e ^' p> a wk {x + 1, x + 2) + 

e(j)e(-u)q- <e ^'? > a ujk (x,x + 2) + £ e(j)e(l)q- <e ^>b lk (x, x + 2), (3.43) 

I 

Eb jj (x,x + 2) = 

e{j)e{uj)q- <e ^ p> a uj {x + l,x + 2) + e{j)e{-uj)q~ <e ^ u ' p> a^{x, x + 2) + 



(-])e(-uj)q <e > +e ^ p> a ju (x, x + 1) + e(-j)e(c;)g <ej "-' p> a^(a;, x + 2) + 



£e(iM0cr <ej+ei ' P> M^ + 2) + £e(-j>(/)g <e ^ ei ' p> M^ + 1)- (3-44) 
i I 

They can be solved modifying the parametrization for nearest neighbors in the 

usual way. The result is: 



0^(1) = r{C l )i2+F aku , b jk {2) = e{j)q-<^>e 2 +r b . k , (3.45) 



a k M) = riC^tfti + Fa^ 6*i(l) = eO')?" <e " P> ^i+^, (3-46) 



where 



^ = -6(^)g-<-^>^,^.=e(fc)g-< e '=^^ fc+ , i o, 

^ = -e(-uj)q <e »> p >Z 2 ,F bjh = e(-k)q< e ^>ti 2 5 k ^ (3.47) 



with ^ = 1,^ 2 £ 2 = 1 resulting from periodic boundary conditions. 
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4 Free boundary conditions 

It is for free boundary conditions, that the Hamiltonian H(Q) commutes with the 
quantum group U q (X n ). As expected, the Bethe states are highest weight states 
of U q (X n ), except some E = states. Since the extension of the BA- procedure 
from the periodic boundary conditions to the free case, follows exactly the lines of 



ref . [ 11 1 , we will only state the results for the sector r = 4. 

Take the block e\ — 62 — 1. The nearest approach constants to be added to the 
now standing waves are the same as in the periodic case, namely equs.( [O0| ),( |3.31|) , 



only jF fc .. (2) is different [II] ) . The BA equations are now: 



where 



C N = n a = 1,2, (4.1) 

b=l,b^a °VSaj?6j 



6) = j- [6 + C 1 - 2A - E ab \ [C 1 + C 1 " 2 A - E ab ] . (4.2) 



and 

E ab = 2Er e " p> + r(gr(c 1 ) + r(6)r(c 1 ). (4.3) 
i 

The only other block is e\ = €2 = 0, E = 0. It is again highly degenerate with 
.F 6ji (2), £1, £ 2 as free parameters^} 

Thus all models have spectra equivalent to the one of the XXZ model. 



5 Conclusion 

We obtained the spectra of quantum spin chain models, arising as representations of 
the Temperley-Lieb algebra associated with quantum groups. The tool is a modified 
version of the coordinate Bethe Ansatz, since the simpler algebraic Bethe Ansatz is 
not immediately available for these models. We find that all models have equivalent 
t Actually there is one more free parameter, called 0:5 in ref. [fllf . 
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spectra, i.e. they differ at most in their degeneracies. The energy eigenvalues are 
given by 

e = f> n (-2 * a + r&jixc 1 )), (5.1) 

n=l 

where — 2A = J2'i q~ <e ' ,p> and the rapidities £ ra are solutions of the BA equations. 

In the sector r we may have p pseudoparticles N w *_i, N u *_i, ■ • • , iV_ w . + i impu- 
rities of the type {u* — 1), (u* — 2), • • • , (— uj* + 1), respectively, such that 

+ 2iV w ._! + ■ • • + (2u* - 1)JV_ W . +1 = r - 2w*p. (5.2) 

Here cj* = u; for A\ and _B n and uu* — u — 1/2 for C n and £)„. 

For example, for periodic boundary conditions^, the total rapidity £ = £1^2 . . . 4p4irop, 4im P 
4p+i4p+2 • • • Cr-a obeys £ N = 1 and the BA equations for E 7^ are: 

where e a b = E^ + 2A and a can be fixed for each allowed case from the equations 
Q5.2p.The Bethe eigenstates are highest weight states of the quantum group, ex- 
cept for the states with energy E = 0, for which this is not always true, in the 
sector r with p impurities and E 7^ 0, the total rapidity £ = £i£ 2 . . . £ p £j m p, Cimp — 



£»+i£p+2 ■ ■ ■ Cr-p obeys £ N = 1 and the BA equations are: 



[Nt2 TT £a (1 + (£a£b) 1 ja e /r a\ 

Sa Simp 11 t (\ \ ( t t \—lt c ' ^ ' 

6=li6 ^ a 4b (i + (4a4feJ 46 - e 
where e = E + 2A. The Bethe eigenstates are highest weight states of the quantum 
group, except for the states with energy E = 0, for which this is not always true. 
These results are expected, but as far as we know, unprooven using only Temperley- 
Lieb algebraic statements as input. 
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